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Abstract 

Consider any supercritical Galton- Watson process which may become extinct with 
positive probability. It is a well-understood and intuitively obvious phenomenon that, 

i-£h on the survival set, the process may be pathwise decomposed into a stochastically 'thin- 

ner' Galton- Watson process, which almost surely survives and which is decorated with 
immigrants, at every time step, initiating independent copies of the original Galton- 
Watson process conditioned to become extinct. The thinner process is known as the 
backbone and characterizes the genealogical lines of descent of prolific individuals in 
the original process. Here, prolific means individuals who have at least one descendant 
in every subsequent generation to their own. 

^^ | Starting with Evans and O'Connell [17], there exists a cluster of literature, |13|. 

|30| [5j [21 [26], describing the analogue of this decomposition (the so-called backbone 
decomposition) for a variety of different classes of superprocesses and continuous-state 

_~ branching processes. Note that the latter family of stochastic processes may be seen 

as the total mass process of superprocesses with non-spatially dependent branching 
mechanism. 

In this article we consolidate the aforementioned collection of results concerning 
backbone decompositions and describe a result for a general class of supercritical su- 
perprocesses with spatially dependent branching mechanisms. Our approach exposes 
the commonality and robustness of many of the existing arguments in the literature. 
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1 Superprocesses and Markov branching processes 

This paper concerns a fundamental decomposition which can be found amongst a general 
family of superprocesses and has, to date, been identified for a number of specific sub-families 
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thereof by a variety of different authors. We therefore start by briefly describing the general 
family of superprocesses that we shall concern ourselves with. The reader is referred to the 
many, and now classical, works of Dynkin for further details of what we present below; see 
for example [6j El El [TQ] . The books of Le Gall [28], Etheridge [H] and Li [29] also serve 
as an excellent point of reference. 

Let E be a domain of ~R d . Following the setting of Fitzsimmons [20], we are interested in 
strong Markov processes, X = {X t : t > 0} which are valued in M.p(E), the space of finite 
measures with support in E. The evolution of X depends on two quantities V and if). Here, 
V = {Vt '■ t > 0} is the semi-group of a diffusion on E (the diffusion killed upon leaving E), 
and ip is a so-called branching mechanism which, by assumption, takes the form 

ip(x,\) = -a(x)X + (3(x)\ 2 + f (e~ Xz -l + \z)n(x,dz), (1) 

J(0,oo) 

where a and /3 > are bounded measurable mappings from E to M and [0, oo) respectively 
and for each x G E, n(x, dz) is a measure concentrated on (0, oo) such that x — > f, Q .(z A 
z 2 )ir(x, dz) is bounded and measurable. For technical reasons, we shall additionally assume 
that the diffusion associated to V satisfies certain conditions. These conditions are lifted 
from Section II. 1.1. (Assumptions 1.1A and LIB) on pages 1218-1219 of [7)H They state 
that V has associated infinitesimal generator 

d 2 d 

L =y a,i j— — - — h y bi— — , 

'--' ' OXiOXj / — ' OXi 
i,j J i 

where the coefficients a^j and bj are space dependent coefficients satisfying: 
(Uniform Elliptically) There exists a constant 7 > such that 
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/] aijUiUj > 7 2J u ~ 

i.l i 



for all x G E and U\, • ■ ■ ua G 



(Holder continuity) The coefficients a^ and 6, are uniformly bounded and Holder contin- 
uous in such way that there exist a constants C > and a G (0, 1] with 

\a,ij(x) - a,ij(y)\, \bi(x) - bi(y)\ < C\x - y\ a 

for all x,y G E. Throughout, we shall refer to X as the (V, "0)-superprocess. 

For each /i G M.p{E) we denote by P M the law of X when issued from initial state X = fi. 
The semi-group of X, which in particular characterizes the laws {P M : \x G Mf(E)}, can be 
described as follows. For each \x G M. F [E) and all / G bp(.E'), the space of non-negative, 
bounded measurable functions on E, 

E M (e- (/ ' Xt) ) = exp j - / u f (x,t)fi(dx) 1 t > 0, (2) 



^^The assumptions on V may in principle be relaxed. The main reason for this imposition here comes in 
the proof of Lemma |6 . 1 1 where a comparison principle is used for diffusions. 



where Uf(x,t) is the unique non-negative solution to the equation 

u f {x,t)=Vt[f](x)- [ ds-V s ty(;u f (;t-s))]{x) xeE,t>0. (3) 

JO 

See for example Theorem 1.1 on pages 1208-1209 of [7j or Proposition 2.3 of [20]. Here we 
have used the standard inner product notation, 

(f,t*) = / f(x)Kdx), 

JE 

for n G M.f{E) and any / such that the integral makes sense. 

Suppose that we define £ = {(l,X t ) = for some t > 0}, the event of extinction. For 
each x G E write 

w(x) = -\og¥ Sx (£). (4) 

It follows from (El) that 



E /1 (e-"< 1 ' x *>) = escp j- / u e (x, t)fi(dx)\ t > 0, 



(5) 



Note that ug(t,x) is increasing in 9 and that P /x ((l,X t ) = 0) is monotone increasing. Using 
these facts and letting — > oo, then i — > oo, we get that 

P M (£) = lim P M ((l,X t ) = 0) = exp (-■ / lim lim Mfl (x, t)//(dx) 1 . (6) 

By choosing fi = 5 X , with x £ E, we see that 

P„(£) = exp (- I w{x)fi{dx)\ . (7) 
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For the special case that ip does not depend on x and P is conservative, (l,X t ) is a 
continuous state branching process. If ^(A) satisfy the following condition: 

dA < oo, 



V>(A) 



then P M almost surely we have £ = {limt_ +00 (l, X 4 ) = 0}, that is to say the event of extinction 
is equivalent to the event of extinguishing, see [2] and [26] for examples. 

By first conditioning the event £ on F t := a{X s : s < t}, we find that for all t > 0, 

The function w will play an important role in the forthcoming analysis and henceforth we 
shall assume that it respects the following property. 

(A): w is locally bounded away from and oo. 



The pathwise evolution of sup erpro cesses is somewhat difficult to visualise on account 
of their realisations at each fixed time being sampled from the space of finite measures. 
However a related class of stochastic processes which exhibit similar mathematical properties 
to superprocesses and whose paths are much easier to visualise is that of Markov branching 
processes. A Markov branching process Z = {Z t : t > 0} takes values in the space M. a (E) of 
finite atomic measures in E taking the form J^" =1 5 Xi , where n e NU {0} and X\, • • • , x n 6 E. 
To describe its evolution we need to specify two quantities, (V,F), where, as before, V is 
the semi-group of a diffusion on E and F is the so-called branching generator which takes 
the form 

F(x,s) = q(x)Y,Pn(x)(s n -s), xeE,se [0,1], (8) 

n>0 

where q is a bounded measurable mapping from E to [0, oo) and, the measurable sequences 
{p n (x) : n > 0}, x E E, are probability distributions. For each v e Ai a (E), we denote by P u 
the law of Z when issued from initial state Z = v. The probability P u can be constructed in 
a pathwise sense as follows. From each point in the support of v we issue an independent copy 
of the diffusion with semi-group V. Independently of one another, for (x,t) gBx [0,oo), 
each of these particles will be killed at rate q(x)dt to be replaced at their space-time point 
of death by n > particles with probability p n (x). Relative to their point of creation, new 
particles behave independently to one another, as well as to existing particles, and undergo 
the same life cycle in law as their parents. 

By conditioning on the first split time in the above description of a (V, F)-Markov branch- 
ing process, it is also possible to show that 

E u (e~ {f > Zt) ) = exp | - / v f (x, t)u(dx) i t > 0, 

where Vf(x,t) solves 

e -v f {x,t) = p[ e -/]( x ) + f ds .p s [F(-,e-^ ( -'*- s) )](x) xeE,t>0. (9) 

Jo 

Moreover, it is known, cf. Theorem 1.1 on pages 1208-1209 of [7], that the solution to 
this equation is unique. This shows a similar characterisation of the semi-groups of Markov 
branching processes to those of superprocesses. 

The close similarities between the two processes become clearer when one takes account 
of the fact that the existence of superprocesses can be justified through a high density scaling 
procedure of Markov branching processes. Roughly speaking, for a fixed triplet, /i, V, ip, one 
may construct a sequence of Markov branching processes, say {Z^ : n > 1}, such that the 
n-th element of the sequence is issued with an initial configuration of points which is taken 
to be an independent Poisson random measure with intensity nfx and branching generator 
F n satisfying 

F n (x,s) = -[ip(x,n(l - s)) + a(x)n(l - s)], xeE.se [0,1]. 

It is not immediately obvious that the right-hand side above conforms to the required struc- 
ture of branching generators as stipulated in OH]), however this can be shown to be the case; 



see for example the discussion on p. 93 of [29]. It is now a straightforward exercise to show 
that for all / G bp(-B) and t > the law of (f,n~ l Z\ n ) converges weakly to the law of 
(/, X t ), where the measure X t satisfies ([2]). A little more work shows the convergence of the 
sequence of processes {n~ 1 Z^ : n > 1} in an appropriate sense to a (V, ^-superprocess 
issued from an initial state /i. 

Rather than going into the details of this scaling limit, we focus instead in this paper 
on another connection between superprocesses and branching processes which explains their 
many similarities without the need to refer to a scaling limit. The basic idea is that, under 
suitable assumptions, for a given (V, ^)-superprocess, there exists a related Markov branch- 
ing process, Z, with computable characteristics such that at each fixed t > 0, the law of Z t 
may be coupled to the law of X t in such a way that, given X t , Z t has the law of a Poisson 
random measure with intensity w(x)X t (dx), where w is given by (j4]). The study of so-called 
backbone decompositions pertains to how the aforementioned Poisson embedding may be 
implemented in a pathwise sense at the level of processes. 

The remainder of this paper is structured as follows. In the next section we briefly 
review the sense and settings in which backbone decompositions have been previously stud- 
ied. Section [3] looks at some preliminary results needed to address the general backbone 
decomposition that we deal with in Sections [U and [HI 

2 A brief history of backbones 

The basic idea of a backbone decomposition can be traced back to the setting of Galton- 
Watson trees with ideas coming from Harris and Sevast'ymiov; cf Harris [23]. Within any 
supercritical Galton- Watson process with a single initial ancestor for which the probability 
of survival is not equal to or 1, one may identify prolific genealogical lines of descent 
on the event of survival. That is to say, infinite sequences of descendants which have the 
property that every individual has at least one descendant in every subsequent generation 
beyond its own. Together, these prolific genealogical lines of descent make a Galton- Watson 
tree which is thinner than the original tree. One may describe the original Galton- Watson 
process in terms of this thinner Galton- Watson process, which we now refer to as a backbone, 
as follows. Let < p < 1 be the probability of survival. Consider a branching process 
which, with probability 1 — p, is an independent copy of the original Galton- Watson process 
conditioned to become extinct and, with probability p, is a copy of the backbone process, 
having the additional feature that every individual in the backbone process immigrates an 
additional random number of offspring, each of which initiate independent copies of the 
original Galton- Watson process conditioned to become extinct. With an appropriate choice 
of immigration numbers, the resulting object has the same law as the original Galton- Watson 
process. 

In Evans and O'Connell p2], and later in Englander and Pinsky [13], a new decompo- 
sition of a supercritical superprocess with quadratic branching mechanism was introduced 
in which one may write the distribution of the superprocess at time t > as the result of 
summing two independent processes together. The first is a copy of the original process con- 
ditioned on extinction. The second process is understood as the superposition of mass from 
independent copies of the original process conditioned on extinction which have immigrated 



'continuously' along the path of an auxiliary dyadic branching particle diffusion which starts 
with a random number of initial ancestors whose cardinality and spatial position is governed 
by an independent Poisson point process. The embedded branching particle system is known 
as the backbone (as opposed to the spine or immortal particle which appears in another re- 
lated decomposition, introduced in Roelly-Coppoletta and Rouault [31] and Evans [IE]). In 
both [17] and [13] the decomposition is seen through the semi-group evolution equations 
which drive the process semi-group. However no pathwise construction is offered. 

A pathwise backbone decomposition appears in Salisbury and Verzani [30], who consider 
the case of conditioning a super-Brownian motion as it exits a given domain such that the 
exit measure contains at least n pre-specified points in its support. There it was found 
that the conditioned process has the same law as the superposition of mass that immigrates 
in a Poissonian way along the spatial path of a branching particle motion which exits the 
domain with precisely n particles at the pre-specified points. Another pathwise backbone 
decomposition for branching particle systems is given in Etheridge and Williams [15] , which 
is used in combination with a limiting procedure to prove another version of Evan's immortal 
particle picture. 

In Duquesne and Winkel [5] a version of the Evans-O'Connell backbone decomposition 
was established for more general branching mechanisms, albeit without taking account of spa- 
tial motion. In their case, quadratic branching is replaced by a general branching mechanism 
ip which is the Laplace exponent of a spectrally positive Levy process and which satisfies the 
conditions < — ?//(0+) < oo and f°° l/^>(£)d£ < oo. Moreover, the decomposition is offered 
in the pathwise sense and is described through the growth of genealogical trees embedded 
within the underling continuous state branching process. The backbone is a continuous-time 
Galton Watson process and the general nature of the branching mechanism induces three 
different kinds of immigration. Firstly there is continuous immigration which is described 
by a Poisson point process of independent processes along the trajectory of the backbone 
where the rate of immigration is given by a so-called excursion measure which assigns zero 
initial mass, and finite life length of the immigrating processes. A second Poisson point 
process along the backbone describes the immigration of independent processes where the 
rate of immigration is given by the law of the original process conditioned on extinguishing 
and with a positive initial volume of mass randomised by an infinite measure. This accounts 
for so-called discontinuous immigration. Finally, at the times of branching of the backbone, 
independent copies of the original process conditioned on extinguishing are immigrated with 
randomly distributed initial mass which depends on the number of offspring at the branch 
point. The last two forms of immigration do not occur when the branching mechanism is 
purely quadratic. 

Concurrently to the work of [5] and within the class of branching mechanisms correspond- 
ing to spectrally positive Levy processes with paths of unbounded variation (also allowing 
for the case that — ip'(0+) = oo), Bertoin et al. [3] identify the aforementioned backbone as 
characterizing prolific genealogies within the underling continuous state branching process. 

Berestycki et al. [2] extend the results of [T7] and [5], showing that for sup erpro cesses with 
relatively general motion and non-spatial branching mechanism corresponding to spectrally 
positive Levy processes with finite mean, a pathwise backbone decomposition arises. The role 
of the backbone is played by a branching particle diffusion with the same motion operator 
as the superprocesses and, like Salisbury and Verzani [30], additional mass immigrates along 
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the trajectory of the backbone in a Poissonian way. Finally Kyprianou and Ren [26] look 
at the case of a continuous-state branching process with immigration for which a similar 
backbone decomposition to [2] can be shown. 

As alluded to in the abstract, our objective in this article is to provide a general backbone 
decomposition which overlaps with many of the cases listed above and, in particular, exposes 
the general effect on the backbone of spatially dependent branching. It is also our intention 
to demonstrate the robustness of some of the arguments that have been used in earlier work 
on backbone decompositions. Specifically we are referring to the original manipulations 
associated with the semi-group equations given in Evans and O'Connell [17] and Englander 
and Pinsky [13], as well as the use the Dynkin-Kuznetsov excursion measure, as found in 
Salisbury and Verzani [30], Berestyki et al. [2] and Kyprianou and Ren [26], to describe the 
rate of immigration along the backbone. 

3 Preliminaries 

Before stating and proving the backbone decomposition, it will first be necessary to describe 
a number of mathematical structures which will play an important role. 

3.1 Localisation 

Suppose that the stochastic process £ = {£ t '■ t > 0} on EU{]}, where f is its cemetery state, 
is the diffusion in E corresponding to the semi-group V. We shall denote its probabilities 
by {U x : x G E}. In the next definition, we shall take bp(E x [0,t]) to be the space of 
non-negative, bounded measurable functions on E x [0,£] with the additional property that 
the value of f(x,s) on E x [0, t] is independent of s, and it is implicitly understood that 
for all functions / G bp(E x [0,t]), we extend their spatial domain to include {f} and set 

/({t,«}) = o. 

Definition 3.1 For any open, bounded set D CC E, and t > 0, there exists a random 
measure X® supported on the boundary of D x [0, t) such that, for all / G bp(_E x [0,t]) 
and \i G AIf(-D), the space of finite measures on D, 

- logE^ (e-<-^>) = f uf(x,t)(t(dx), t > 0, (10) 

where u?(x,t) is the unique non-negative solution to the integral equation 



u 



f(x,t) = U x [f(Ur^tAr D )]-U x 



tAT 1 



^(&,«/(&,*-s))ds 



'ir 



and r D = inf{£ > 0,£ f G D c }. Note that, here, we use the obvious notation that (f,Xf) = 
Id(Dx\ot)) fiXy^-X-PidXyds). Moreover, with a slight abuse of notation, since their effective 
spatial domain is restricted to D U {+} in the above equation, we treat ip and u? as functions 
in bp(E x [0, £]) and accordingly it is clear how to handle a spatial argument equal to t, as 
before. In the language of Dynkin [9], X® is called an exit measure. 
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Now we define a random measure X® on D such that (f,Xf) = (f,Xf) for any / G 
bp(-D), the space of non-negative, bounded measurable functions on D, where, henceforth, 
as is appropriate, we regard / as a function defined on E x [0, oo) in the sense that 



f(x), x e D 

0, xeE\D. 



Then for any / G bp(D) and pL G M f (D), 

-logE M (e^ x ">) = j uf(x,t)fi(dx) 



(12) 



t > 0, 



(13) 



where u?(x,t) is the unique non-negative solution to the integral equation 



uf(x,t) = n x [f(£ t y,t<T D ]-ii x 



tAT L 



1>{Zs,u?{Z a ,t-s))ds 



x ED. (14) 



As a process in time, X D = {X^ : t > 0} is a superprocess with branching mechanism 
4>(x, A)lu(x), but whose associated semi-group is replaced by that of the process £ absorbed 
on 5.D. Similarly, as a process in time, X D = {X® : t > 0} is a superprocess with branching 
mechanism -0(x, A)1d(x), but whose associated semi-group is replaced by that of the process 
£ killed upon leaving D. One may think of Xf as describing the mass at time t in X which 
historically avoids exiting the domain D. Note moreover that for any two open bounded 
domains, D\ CC D 2 CC E, the processes X Dl and X° 2 (and hence X Dl and X D2 ) are 
consistent in the sense that 



xr = w 



(15) 



for all t > (and similarly X? 1 = (X t D2 ) Dl for all t > 0). 



3.2 Conditioning on extinction 

In the spirit of the relationship between ( ITOi) and ( ITTj) . we have that u> is the unique solution 
to 

^(6»,«>(6))ds 



w(x) = n x [M;(£ iAr D)] - n a 



x e D. 



(16) 



for all open domains D GG E. Again, with a slight abuse of notation, we treat w with its 
spatial domain E U {f} as a function in bp(.E x [0,t]) and w(f) := 0. From Lemma 1.5 in 
[7] we may transform ( Tl6|) to the equation 



wfx) = IL, 



W (^Ar D )eXp<^ - / 



WsMZs)) 



'0 ^fe 

which shows that for all open bounded domains D, 



ds 



x G D, 



w(£t A r«)eXp 



t/\T L 



W{£s) 



ds 



t > 0, 



;i7) 



is a martingale. 

The function w can be used to locally describe the law of the superprocess when condi- 
tioned on global extinction (as opposed to extinction on the sub-domain D). The following 
lemma outlines standard theory. 

Lemma 3.2 Suppose that fi £ Aip(E) satisfies (w, /i) < oo (so, for example, it suffices that 
fj, is compactly supported) . Define 

p;(-) = p^(-|£). 

Then for any f £ bp(E x [0, £]) wi& the additional property that the value of f(x,s) on 
E x [0,i] is independent of s and \x £ M. F {D), 

-logE; ( e -<'^>) = | uf *(M)M<H 

where Ur'*(x,t) = uf +w (x,t) — w(x) and it is the unique solution of 



n 



-D,* 



»'*(x,t)=Il x [f(Z tAT n)]-Il a 



tAr D 



x £ D, 



ip*(Z, s ,u f '*(Z, sl t-s))ds 

where ip*(x, A) = ip(x, A + w(x)) — ip(x,w(x)), restricted to D, is a branching mechanism of 
the kind described in the introduction and for each \i £ Aip(E), (X, P*) is a superprocess. 
Specifically, on E, 

iP*(x,X) = -a*(x)\ + f3(x)\ 2 + f (e~ Xz -l + \z)7i*(x,dz), (19) 

J{0,oo) 

where 

a*{x) = a{x) - 2f3(x)w{x) - f (1 - e- w(x)z )zn(x, dz) 

J(0,oo) 

and 

tt*(x, dz) = e- w{x)z ir(x, dz) on E x (0, oo). 

Proof: For all / £ bp(d(D x [0,t))) we have 

E;(e-^) = E IM (e-^\£) 

= e^%(e-^ D h £ ) 

= e^E^e-V'^Exnile)) 



N 



= e^%(e~V +w ' x ^) 

_ e -{uf +w (-,t)-w,tJ,) _ 

Using ffTTT) and ffT6l) then it is straightforward to check that Uf '*(x, t) = uf +w (x, t) — w(x) is 
a non-negative solution to (|T8|) . which is necessarily unique. The proof is complete as soon 
as we can show that tp*(x, A), restricted to D, is a branching mechanism which falls into 
the appropriate class. One easily verifies the formula (fT9"|) and that the new parameters a* 
and 7r*, restricted to D, respect the properties stipulated in the definition of a branching 
mechanism in the introduction. □ 
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Corollary 3.3 For any bounded open domain D CC E, any function f G bp(D) and any 
ji G A4f(D) satisfying (w,fj) < oo, 



logE;(e~ 



i.f,x t D ) 



D,* 



D 



Ur '*(x, t)/i(dx) 



D,*i 



where Uj '*(x,t) = Uj +W (x,t) — w(x) and it is the unique solution of 

u J /'*(x,t) = U x [f(^);t<r D ]-U x 
where ip* is defined by (TH)]) . 



^*(6,«/'*(6,*-s))ds 



x G D, 



(20) 



3.3 Excursion measure 

Associated to the law of the processes X, are the measures {N x : x G E}, defined on the 
same measurable space as the probabilities {P| : x G E} are defined on, and which satisfy 



N:.(l - e-<^>) = -logE^(e-<^>) = u* f (x,t), 



(21) 



for all / G bp(-E) and t > 0. Intuitively speaking, the branching property implies that P| is 
an infinitely divisible measure on the path space of X, that is to say the space of measure- 
valued cadlag functions, ID>([0, oo) x Ai(E)), and ( |2T|) is a 'Levy-Khinchine' formula in which 
N* plays the role of its 'Levy measure'. Such measures are formally defined and explored in 
detail in [12]. 

Note that, by the monotonicity property, for any two open bounded domains, Di CC 
D 2 CC E, 

(f,X^)<{f,X^) N;-a.e., 

for all / G bp(-Di) understood in the sense of (fl2l) . x G D\ and t > 0. Moreover, for an open 
bounded domain D and / as before, it is also clear that N*(l — e~^' x± ') = u, '*(x, t). 

The measures {N* : x G E} will play a crucial role in the forthcoming analysis in order 
to describe the 'rate' of a Poisson point process of immigration. 



3.4 A Markov branching process 

In this section we introduce a particular Markov branching process which is built from 
the components of the (ip, P)-superprocess and which plays a central role in the backbone 
decomposition. 

Recall that we abuse our notation and extend the domain of w with the implicit under- 
standing that w(]) = 0. Note, moreover, that thanks to (IT7|) . we have that, for x G E, 

w(x)~ 1 w(^t) exp < — J ip(£ s , w(£ s ))/w($ iS )ds > is in general a positive local martingale and 

hence a supermartingale. For each t > 0, let T{ = a(£ s : s <t). Let ( = inf{£ > : £ 4 G {f}} 
be the life time of £. The formula 



dig 

dIL 



^ 



wfo 



w[x) 



cxp 



V'fe^fe)) 



w(6 



d.s 



on{t<C}, t>0,xe£, 



f22) 
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uniquely determines a family of (sub-)probability measures {n™ : x G E}. It is known that 
under these new probabilities, £ is a right Markov process on E; see [32], Section 62. We 
will denote by V w the semi-group of the E U {f}- valued process £ whose probabilities are 
{IT- : x e E}. 

Remark 3.4 The equation ( TT6|) may formally be associated with the equation Lw(x) — 
ip(x,w(x)) = on E, and the semi-group V w corresponds to the diffusion with generator 

L™ := L w - w~ l Lw = L W - w"V(", w), 

where L w u = w~ 1 L(wu) for any u in the domain of L. Intuitively speaking, this means that 
the dynamics associated to V w , encourages the motion of £ to visit domains where the global 
survival rate is high and discourages it from visiting domains where the global survival rate 
is low. (Recall from © that larger values of w(x) make extinction of the (ip, P)-superprocess 
less likely under Fg x .) 

Henceforth the process Z = {Z t : t > 0} will denote the Markov branching process whose 
particles move with associated semi-group V w . Moreover, the branching generator is given 
by 

E(x,s) = q(x)Y,Pn(x)(s n -s), (23) 

n>0 

where 

/ \ ,// / xx ib(x,w(x)) /nA . 

q(x) = ij/(x, w(x)) - ^ ' ; JJ , 24 

w{x) 

Po(x) = pi(x) = and for n > 2, 



Here we use the notation 

x e E. 



d 
ip'(x,w(x)) := —ip(x,\) 



\=w(x) 



Note that the choice of q(x) ensures that {p n (x) '■ n > 0} is a probability mass function. In 
order to see that q(x) > for all x 6 E (but q ^ 0), write 



q(x) = P(x)w(x) + —r^ / (l-e- w{x)z {l + w{x)z))<ir{x,dz) (25) 

w{x) y (0jOo) 

and note that (3 > 0, w > and 1 — e~ Az (l + Xz), A > 0, are all non-negative. 

Definition 3.5 In the sequel we shall refer to Z as the (V w , F)-backbone. Moreover, in the 
spirit of Definition 13.11 for all bounded domains D and t > 0, we shall also define Zf to be 
the atomic measure, supported on d(D x [0, £)), describing particles in Z which are first in 
their genealogical line of descent to exit the domain D x [0, £). 
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Just as with the case of exit measures for superprocesses, we define the random measure, 
Z D = {Zf : t > 0}, on D such that (/, Zf) = (/, Zf) for any / G bp(-D), where we remind 
the reader that we regard / as a function defined on E x [0, oo) as in 012|) . As a process in 
time, Z D is a Markov branching process, with branching generator which is the same as in 
(|23|) except that the branching rate q{x) is replaced by q D (x) := q(x)ln(x), and associated 
motion semi-group given by that of the process £ killed upon leaving D. Similarly to the 
case of superprocesses, for any two open bounded domains, D\ CC D 2 CC E, the processes 
Z Dl and Z D2 (and hence Z Dl and Z D2 ) are consistent in the sense that 

z* = (zf 2 ) Dl 

for all t > (and similarly Zf 1 = (Zf 2 ) Dl for all t > 0). 

4 Local backbone decomposition 

We are interested in immigrating (V, ^*)-superprocesses onto the path of an (V w , F)-backbone 
within the confines of an open, bounded domain D CC E and initial configuration v G 
M. a (D), the space of finite atomic measures in D of the form X^=i ^' wnere n6lU {0} 
and Xi, • • • ,x n G D. There will be three types of immigration: continuous, discontinuous 
and branch-point immigration which we now describe in detail. In doing so, we shall need 
to refer to individuals in the process Z for which we shall use classical Ulam-Harris nota- 
tion, see for example p290 of Harris and Hardy [22]. Although the Ulam-Harris labelling 
of individuals is rich enough to encode genealogical order, the only feature we really need 
of the Ulam-Harris notation is that individuals are uniquely identifiable amongst T, the set 
labels of individuals realised in Z. For each individual u G T we shall write b u and d u for 
its birth and death times respectively, {z u (r) : r G [fe u ,rf u ]} for its spatial trajectory and N u 
for the number of offspring it has at time d u . We shall also write T D for the set of labels of 
individuals realised in Z D . For each u G T D we shall also define 

r^ = inf{s G [b u , d u ], z u (s) G D c }, 

with the usual convention that inf := oo. 

Definition 4.1 For v G M. a (D) and \x G M.p{D), let Z D be a Markov branching process 
with initial configuration u, branching generator which is the same as in (l2"5|) . except that the 
branching rate q(x) is replaced by q D (x) := q(x)lr>(x), and associated motion semi-group 
given by that of V w killed upon leaving D. Let X D) * be an independent copy of X D under 
P* Then we define the measure valued stochastic process A D = {Af : t > 0} such that, for 

a? = x? >* + ir* + ir + 1?'\ (26) 

where I D ' W = {/f' N * : t > 0}, I D ' r = {if r : t > 0} and I D " = {lf' v : t > 0} are defined 
as follows. 

i) (Continuum immigration:) The process 1°'^ is measure-valued on D such that 

rD,N* _ V"^ V^ v-(D,l,u,r) 

ueT D b u <r<tAd u Ar^ 
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where, given Z D , independently for each u G T D such that b u < t, the processes 
X( D ' 1 ' u ' r ) are independent copies of the canonical process X D , immigrated along the 
space-time trajectory {(z u (r),r) : r G (b u ,t A d u A t®]} with rate 

dr x 2 j 9(z B (r))dJC B(r) . 
ii) (Discontinuous immigration:) The process I D ' ¥ is measure- valued on D such that 

rD,P* __ V^ V^ v-(D,2,u,r) 

uET D b u <r<tAd u ArO 

where, given Z D , independently for each u G T D such that b u < t, the processes 
X( D ' 2 ' u ' r ) are independent copies of the canonical process X D , immigrated along the 
space-time trajectory {(z u (r),r) : r G (b u ,t A d u A r^]} with rate 

dr x / ye-^-MJM^r), dy) x dP^ 

</j/G(0,oo) 

iii) (Branch point biased immigration:) The process I D,ri is measure-valued on D such 

that 

T D, V _ V^ -, v(-D,3,ti) 

wGT D 

where, given Z D , independently for each w G 7"° such that d u < t A rf, the processes 
j^(-D,3,n) are independent copies of the canonical process X D issued at time d u with law 
Py 5 such that, given u has n > 2 offspring, the independent random variable F M 

has distribution r) n (z u (r) , dy) , where 

T] n (x,dy) = -—-. - \p(x)w 2 (x)5 (dy)l {n=2} + w(x) n ^-e- w ^yTT(x,dy) 

q{x)w{x)p n {x) { n\ 

'(27) 

It is not difficult to see that A D is consistent in the domain D in the sense of (TT5|) . 
Accordingly we denote by P( /i ,y) the law induced by {A®,D G 0(E),t > 0}, where O(E) is 
the collection of bounded open sets in E. 

The so-called backbone decomposition of (X D ,F^) for ji G M. F [D) entails looking at the 
process A D in the special case that we randomise the law Pq^) by replacing the deterministic 
choice of v with a Poisson random measure having intensity measure w(x)fi(dx). We denote 
the resulting law by P^. 

Theorem 4.2 For any fi G Aip(D), the process (A^P^) is Markovian and has the same 
law as (X D ,P^). 
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5 Proof of Theorem 14.2 



The proof involves several intermediary results in the spirit of the non-spatially dependent 
case of Berestycki et al. [2]. Throughout we take D as an open bounded domain in E. Any 
function / defined on D will be extend to E by defining / = on E \ D. 

Lemma 5.1 Suppose that [i E M. F (D), v E A4 a (D), t > and f E bp(-D), We have 
E M (e-tt^+tf^KZf :s<t}) = exp |- j\cf>(., uf *(•, t - *)), Zf )ds} , 



where 

<f>(x,\) = 2(3{x)\ + 

'(0,oo) 



x,X) = 2/3(x)X + {l-e~ Xz )zir(x,dz), xED,\>0. (28) 

J(0,oo) 



Proof: We write 

</, ir* + iD = E E </. ^"' r) > + E E </. *£ 2, " r) >- 

Hence conditioning on Z D , appealing to the independence of the immigration processes 
together with Campbell's formula and that N*(l — e~^' Xs •*) = u , '*(x, s), we have 

E ( ^)(e-^ D ' N *)|{Zf :s<t}) 

t/\d u ATf> 



ueT D 
On the other hand 

E M (e-^ F \ Z n :s < t}) 

-tAduArJ? poo 



{ptAa u AT~ 
" E 2 / P(z u (r)).N* Zu{r) (l-e^ x ^)dr 

u£T D K 

frtAduAT° "] 

- E 2 J P(z u (r))uf>*(z u (r),t-r)dr\. (29) 



iptACt u AT~ poo 
E / / ye- w{Mr))y n(z u (r),dy)E; Szu{r) (l-e-^ x ^)dr 

ueTD Jb u JO 

{ptAduAT® poo ~\ 

"E / / ^-e^ u f'* {Zu{r) ' t ~ r)y )ye- wMr))y 7i{z u {r),dy)dr\. (30) 

ueTD Jb u Jo J 

Combining (|2"9"|) and (13"Uj) the desired result follows. D 



Lemma 5.2 Suppose that the real-valued function J(s, x, A) defined on [0,T)xDxK satisfies 
that for any c > there is a constant A(c) such that 

| J(s, a;, Ai) - J{s, x, A 2 )| < A(c)|Ai - A 2 |, 
14 



for all s G [0,T), x G D and Ai,A 2 G [— c, c]. T/ien for any bounded measurable function 
g(s,x) on [0,T) x D, the integral equation 

v(t, x) = g(t, x)+ f U x [J(t - s, &, v(t - s, £ s )); s < r D ] dx, t G [0, T), 
Jo 

/ms ai mos£ one bounded solution. 

Proof: Suppose that v\ and v% are two solutions, then there is a constant c > such that 
—c<Vi, V2 < c and 

\\vi-v 2 \\{t) <A{c) [ \\v 1 -v 2 \\{s)ds, 
Jo 

where \\v\ —V2\\(t) = sup xeD \vx(t,x) — V2(t,x)\, t G (0,T). It follows from Gronwall's lemma 
(see, for example, Lemma 1.1 on page 1208 of [7j) that ||i>i — t>2||(£) = 0,t G [0, T). 

Lemma 5.3 Fix t > 0. Suppose that f,h G bp(Z?) and (? s (x) is jointly measurable in 
(x, s) £ D x [0,t] and bounded on finite time horizons of s such that g s (x) = for x G D c . 
Then for any \i G J\A. F [D), x G D and t > 0, 



e -W{x,t) . = E 



(lJ.,6x) 



exp 



{^.,z?)a8-(f,i?«)-(h,zi>) 



where e H/ ( x ' t ) is the unique [0, l]-valued solution to the integral equation 



w(x)e 



-W(x,t) 



Tl x [w(Z t )e- h ^\t<T D ] 



+IL 



t/\T D 



[H t - S (£s, -w(£ 8 )e 



-Wfo,t-s)> 



- w(£a)e 



-Wfa,t-s) 



gt-s(Cs 



-^ s ,w(Q)e- w ^^]d± 



(31) 



for x G D, where 

POO 

H t „ s {x, A) = q(x)X + (5{x)\ 2 + / (e~ Xy - 1 + Aw) e - (w(x)+u / '"^^Kix, dy), x G £>, 

Jo 

and q{x) was defined in (24\ ). 

Proof: Following Evans and O'Connell [T7] it suffices to prove the result in the case when 
g is time invariant. To this end, let us start by defining the semi-group V h,D by 



v^[k](x) = nJe-C h ^k(t tAT ») 



n^ e -/oM6)ds fc(6) . t<r Dj D 



where, for convenience, we shall write 



v t D ik]=-prik] 



(32) 



(33) 
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Define the function x( x ) = ip(x,w(x))/w(x). Conditioning on the first splitting time 
in the process Z D and recalling that the branching occurs at the spatial rate q D (x) = 
l^(a;)('i/' / (x, w(x)) — x( x )) we get that for any x G D, 

e- w ^ l) -^V? +q+x > D [we~ h }{x) 



+n: 



t/ST 



W{X) 

D 



(34) 



exp I - / (g + q)(£ r )dr 

q(L) Y,P^s)e- nW ^- s) [ r,n(ts, dy)e~^?'^ t -^ 1 & 

n>2 J (°>°°) J 

From (1271) we quickly find that for x G D 

y> n (x)e-^'*- s ) / Vn (x,dy)e- yu f'*^ t - s) 

_^n J(0,oo) 

{H t „ s (x, -w(x)e- w{x ' t - s) ) + w(x)q(x)e- w{x ' t - s) } . 



(35) 



n>2 



q(x)w(x 
Using the above expression in ( )34l) we have that 



w(x)e 



-W(x,t) 



+n a 



v\ 



g+q+x,D 



we h ](x) 



/ exp I- (g + q + x)(£r)dr\ 



[(Ht-fa-wiQe 



-w(e«,t-s)> 



«' 



(e s )g(e s )e- W(es '*- s) )l ds 



Now appealing to Lemma 1.2 in Dynkin [10] and recalling that %(•) — ' l P('> w ('))/ w (') on -D> 
we may deduce that for any igfl, 



■w 



(x)e~ w{x ' t] = V t D lwe~ h ](x) 



+n a 



tAr 1 



[#t_ s (£ s ,-u;(£ s )e 



-W(£,,t-s)^ 



^(6)^(6)e 



-W(6,t-s) 



-^fe,^fe))e- W(6 ' 4 - s) ]d S 



(36) 



as required. Note that in the above computations we have implicitly used that w is uniformly 
bounded away from and ooonD. 

To complete the proof we need to show uniqueness of solutions to fl36l) . Lemma [531 offers 
sufficient conditions for uniqueness of solutions to a general family of integral equations which 
includes fl36|) . In order to check these sufficient conditions, let us define w D = sup yeD w(y). 
Thanks to Assumption (A) we have that < w D < oo. For s > 0, x G D and A G [0,uJ D ], 
define the function J(s,x, A) := [H s (x, —A) — (g(x) + x(x))X]. We rewrite ( 1361) as 



w(x)e 



V?[we- h ](x)+ f U x [J(t-s^ s ,w^s)^ Wi$s ^ s) y,s<r D ]ds. 
Jo 



-W(x,t) _ vDr-h 
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Lemma [5^2l tells us that (1361) has a unique solution as soon as we can show that J is continuous 
in s and that for each fixed T > 0, there exists a K > (which may depend on D and T) 
such that 

supsup \J(s,y,Xi) - J(s,y,X 2 )\ < K\\ ± - A 2 |, Ai, A 2 G (0,w D ]. 

s<T y£D 

Recall that g(y) is assumed to be bounded, moreover, Assumption (A) together with the 
fact that 

sup\\a(y)\+P(y)+ (z A z 2 )<ir(y,dz) 1 < oo (37) 

»/GD L J(0,oo) J 

also implies that x is bounded on D. Appealing to the triangle inequality, it now suffices to 
check that for each fixed T > 0, there exists a K > such that 

snpsnp\H s (y,-\ 1 )-H s (y,-X 2 )\<K\\ 1 -\ 2 \, \ u A 2 G (0, «J D ]. (38) 

s<TyeD 

First note from Proposition 2.3 of Fitzsimmons [20] that 

sup sup M f '*(x, s) < oo. (39) 

Straightforward differentiation of the function H s {x, —A) in the variable A yields 



H.(x, -A) = q(x) - 2(3(x)X + [ (1 - e Az )e" (u;(a;)+n ?'* (x ' s)) ^7r(x, dz). 



d_ 
9X j (o,oo 

Appealing to ( l37|) and ( !39l) it is not difficult to show that the derivative above is uniformly 
bounded in absolute value for s < T, x G D and A G [0,w D ], from which ( )38l) follows by 
straightforward linearisation. The proof is now complete. □ 



Theorem 5.4 For every /i G M. F (D), v G A4 a (-D) and f,h G bp(-D) 

E^j (e-tfA^-W)) = e -<«f* (-,*).m>-(^ (-A-), (40 ) 



where e~ v f> h ^ x '' is the unique [0,l]-solution to the integral equation 
)e -vf, h (x,t) = n x [w(tt)e- h M;t<T D ] 



w(x) 



+n x 



tAr D 



/ 





.(41) 



Proof: Thanks to Lemma [3.21 it suffices to prove that 

where I D := J D > N * + / D . p * _|_ I D >V } anc [ v d so i ves (J4JJ). Putting Lemma 15.11 and Lemma 
15.31 together we only need to show that, when gt- s {') — (fii^Uf '*(■>£ — s)) (where <j) is given 
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by (|28|) ). we have that exp{— W(x, t)} is the unique [0, l]-valued solution to f|4TT) . Again 
following the lead of [2], in particular referring to Lemma 5 there, it is easy to see that on D 

J w(-) 

= r(; w{-)e- w <--*-> + «?'*(-, t-s))- p(., uf'*(-, t-s)), 

which implies that exp{— W(x, t)} is the unique [0, l]-valued solution to (J4"T|) . D 

Proof of Theorem 14. 2t . The proof is guided by the calculation in the proof of Theorem 
2 of [2j. We start by addressing the claim that (A D ,P M ) is a Markov process. Given 
the Markov property of the pair (A D ,Z D ), it suffices to show that, given Af , the atomic 
measure Z® is equal in law to a Poisson random measure with intensity w(x)A®. Thanks 
to Campbell's formula for Poisson random measures, this is equivalent to showing that for 
all h E bp(.D), 

which in turn is equivalent to showing that for all f,h 6 bp(-D), 

E At (e-^A D M^?>) = jy e -K(i-^)+/,A?>). (42) 

Note from (HUj) however that when we randomize v so that it has the law of a Poisson random 
measure with intensity w(x)fi(dx), we find the identity 



EJe 



J i' 



-</,a?>-<MP>) = exp (-uf*(.,t) -w(l- e-^(-^),/i) . 



Moreover, if we replace / by w • (1 — e h ) + / and h by in ()40|) and again randomize v so 
that it has the law of a Poisson random measure with intensity w(x)fi(dx) then we get 

E, ( e -<-a-- h H/A D >) = exp (-<:* 1 _ e _ h)+/ (-^) - «, ■ (l - exp {-r£ (1 _ B _ h)+/)0 }) ,/z 

These last two observations indicate that (142|) is equivalent to showing that, for all /, h as 
stipulated above and t > 0, 

«f *(x,t) + w(x)(l - e -^ (x ' 4) ) = i# ( \_ e _ h)+/ (z,i) + w(x)(l - e "^-(i-e-") + /,o ( "' 4) ). (43) 

Note that both left and right-hand side of the equality above are necessarily non-negative 
given that they are Laplace exponents of the left and right-hand sides of (1421) . Making use 
of (fT8l) . (TT6l) . and (14TT) . it is computationally very straightforward to show that both left and 
right-hand side of (|43j) solve (fTTj) with initial condition / + w(l — e~ h ), which is bounded in 
D. Since (1.2) has a unique solution with this initial condition, namely u 1 ? +w , 1 _ e _ h Jx,t), we 
conclude that (j4"51) holds true. The proof of the claimed Markov property is thus complete. 
Having now established the Markov property, the proof is complete as soon as we can 
show that (A D ,P M ) has the same semi-group as (X D ,F^). However, from the previous part 
of the proof we have already established that when f,h£ bp(£)), 



e m e -^n-tf>An = e -^- c -k )+f ^-»»> = e m e -u+^-e-^n . (44) 



In particular, choosing h = we find 



£ > 0, 



which is equivalent to saying that the semi-groups of (A , P M ) and (X , P M ) agree. □ 



6 Global backbone decomposition 

So far we have localized our computations to an open bounded domain D. Our ultimate 
objective is to provide a backbone decomposition on the whole domain E. To this end, let 
D n be a sequence of open bounded domains in E such that D n f E. Let X Dn , A Dn and Z Dn 
be defined as in previous sections with D being replaced by D n . 

Lemma 6.1 For any h, f G bp(E') wift compact support and any // G .M.f(.E'), we aai>e 
£aa£ /or any £ > 0,each element of the pair {((h, Z^ n ), (/, Af")) : s < £} pathwise in- 
creases P ^-almost surely as n — > oo. T7ie limiting pair of processes, here denoted by 
{((h,Zf n ), (/, Af n )) : s < t}, are swca £aa£ (f,Af n ) is equal in law to (f,X t ) and, given 
A™ m ; the law of Z™ m is a Poisson random field with intensity w(x)Af im (dx). Moreover, Z mm 
is a (V W ,F) branching process with branching generator as in [2&\l and associated motion 
semi-group given by that of the process £. 

Proof: Appealing to the stochastic consistency of Z D and A D in the domain D, it is 



D 



clear that both {h, Z t n , ) and (/, A t n ) are almost surely increasing in n. It therefore follows 



that the limit as n — > oo exists for both (h, Z t n ) and (/, A t ™), P M -almost surely. In light 
of the discussion at the end of the proof of Theorem 14.21 the distributional properties of the 
limiting pair are established as soon as we show that 



logE^ [e 



(h,zr n )-(fAf n ) 



u, 



L ,(i_ e -h )+/ (x,£)/i(dx), t > 0. 



(45) 



Assume temporarily that /i has compact support so that there exists an n G N such that 
for n > n we have that supp/i C D n and h = f = on D c n . Thanks to (|44|) and monotone 
convergence (|45[) holds as soon as we can show that u® n "f u g for all g G bp(.E') satisfying 
g = on D c n for n > n . By (TT3|) and (Tl4|) . we know that u^ n (x, t) is the unique non- negative 
solution to the integral equation 



u?«(x,t)=U x \<,(£ tATD )]-IL B 



t/\T Dn 



^ s ,uf n (^,t-s))ds 



(46) 



Using Lemma 1.5 in [7] we can rewrite the above integral equation in the form 

"tAT Dn 



«f"(x,t) 



n. 



g(£tAr D )exp 



afc)ds 



-n. 



*Ar 0n 



expf / a(^)drj [V>(&, <"(&, *-«)) + «(£>£"(&,£ - s))] ds 



(47) 
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Since g = on D c n for n > n , we have 

-tf\m n 



n. 



#(&AT D J eX P 



a(Qds 



n, 



#(&) exp / a(Qds ; t < r Dn 



which is increasing in n. By the comparison principle, uP n is increasing in n (see Theorem 
3.2 in part II of [7]). Put u 9 = lim^^oo u® n . Note that i/)(x, A) + a(x)A is increasing in A. 
Letting n — > oo in (jUJ), by the monotone convergence theorem, 



where 



u g {x,t) = V?g{x) 



V?f = n * 



n* [v. 

Jo 



u g (-, t-s)) + a(-)u g (-, t - s))] ds, 



9(&)exp ( / «fe)ds 



, 9 e bp(E), 



which in turn is equivalent to 

u g (x, t) = V s g{x) - n x . / V s tp{-, u g (; t - s))ds. 

Jo 

Therefore, u g and u g are two solutions of ([3]) and hence by uniqueness they are the same, as 
required. 

To remove the initial assumption that \i is compactly supported, suppose that \i n is 
a sequence of compactly supported measures with mutually disjoint support such that 
fj, = J2k>i^k- By considering (H5l) for Y^k=i^k and taking limits as n | oo we see that 
(|45p holds for /x. Note in particular that the limit on the left hand side of (I45p holds as a 
result of the additive property of the backbone decomposition in the initial state \x. D 



Note that, in the style of the proof given above (appealing to monotonicity and the max- 
imality principle) we can easily show that the processes X Dn '*, n > 1, converge distribution- 
ally at fixed times, and hence in law, to the process (X, P*); that is, a (V, ^*)-superprocess. 
With this in mind, again appealing to the consistency and monotonicity of the local back- 
bone decomposition in the size of domain the following, our main result, follows as a simple 
corollary of Lemma 16.11 

Corollary 6.2 Suppose that // G J\4p(E). Let Z be a (V w ,F)-Markov branching process 
with initial configuration consisting of a Poisson random field of particles in E with intensity 
w(x)fi(dx). Let X* be an independent copy of (X, P*). Then define the measure valued 
stochastic process A = {A t : t > 0} such that, for t > 0, 



A t = x; + i? +/f +/, 



t ■ 



(48) 



where I w = {if : t > 0}, J p * = {If : t > 0} and 1*1 = {i? : t > 0} are defined as follows, 
i) (Continuum immigration:,) The process I w is measure-valued on E such that 



T = £ E 



A" 



(l,u,r) 



t-r 



uGT b u <r<tAd u 

20 



where, given Z , independently for each u G T such that b u < t, the processes X^ l,u ' T ^ 
are independent copies of the canonical process X, immigrated along the space-time 
trajectory {(z u (r),r) : r G (b u ,t A d u ]} with rate 

drx2/3(z u (r))dW Mr) . 
ii) (Discontinuous immigration:^ The process I ¥ is measure-valued on E such that 



u£T b u <r<tAd u 



r) 



where, given Z , independently for each u G T such that b u < t, the processes X^ 2,u,r " > 
are independent copies of the canonical process X, immigrated along the space-time 
trajectory {(z u (r),r) : r G (b u , t A d u }} with rate 

dr x f ye- w ^ r ^7r(z u (r), dy) x dP^ 

Hi) (Branch point biased immigration:^ The process I v is measure-valued on E such 
that 

u<=T D 

where, given Z , independently for each u G T such that d u < t, the processes X^ 3 '^ are 
independent copies of the canonical process X issued at time d u with law ^y u &, such 
that, given u has n > 2 offspring, the independent random variable Y u has distribution 
r) n (z u (r),dy), where r) n (x,dy) is defined by (F7\ ). 

Then (A, P M ) is Markovian and has the same law as (X, P M ). 
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